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Recent Observations of Supernova Remnants with VERITAS 
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Abstract: Supernova remnants (SNRs) are widely considered the most likely source of cosmic rays below the knee 
(10 15 eV). Studies of GeV and TeV gamma-ray emission in the vicinity of SNRs, in combination with multi-wavelength 
observations, can trace and constrain the nature of the charged particle population believed to be accelerated within SNR 
shocks. They may also speak to the diffusion and propagation of these energetic particles and to the nature of the acceler- 
ation mechanisms involved. We report here on recent observations of SNRs with VERITAS, including the discoveries of 
VHE gamma-ray emission from from G120. 1+1.4 (Tycho's SNR) and from the northwest shell of G78.2+2. 1 (gamma-ray 
source VER 12019+407, which was discovered as a consequence of the VERITAS Cygnus region survey). 
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1 Introduction 

Cosmic rays with energies up to 10 15 eV are believed 
to be accelerated at strong shocks of supernova remnants 
(SNRs). Current space-borne (Fermi) and ground-based 
(VERITAS, HESS, MAGIC) gamma-ray telescopes allow 
us to study cosmic -ray accelerators via secondary gamma- 
ray production, either via the inverse Compton process 
(electrons) or via the decay of neutral pions (tt° — > 77) 
produced when high-energy cosmic-ray ions collide inelas- 
tically with the ambient medium. Studies of GeV and TeV 
gamma-ray emission in the vicinity of SNRs, in combi- 
nation with multi-wavelength observations, can trace and 
constrain the nature of the charged particle population be- 
lieved to be accelerated within SNR shocks. They may 
also speak to the diffusion and propagation of these ener- 
getic particles and to the nature of the acceleration mech- 
anisms involved. We discuss here recent detections of two 
TeV gamma-ray sources by VERITAS, one identified with 
Tycho's SNR, the other potentially identified with SNR 
G78.2+2. 1, that have the potential to further illuminate this 
long-standing question. 

2 VERITAS 

The VERITAS instrument is an array of four 12-m imag- 
ing atmospheric Cherenkov telescopes located at the Fred 
Lawrence Whipple Observatory in southern Arizona. Each 
telescope is equipped with a 3.5° field-of-view (FOV). The 
instrument has good (15-25%) energy resolution over an 
energy range of 100 GeV - 30 TeV and excellent angu- 
lar resolution (Rgs < 0.1°), and achieves a sensitivity of 



1% of the Crab Nebula flux in approximately 25 hours. 
For more details on the VERITAS instrument and analy- 
sis techniques, please refer to ITT1IT21 . 

3 Tycho's SNR 

The SNR, G 120. 1+1. 4, also known as Tycho's SNR, is the 
relic of a Type la supernova lfl4l observed in 1572. Tycho's 
SNR is young among Galactic SNRs (438 years). It shows 
a distinct shell-like morphology in radio [5 1 and strong non- 
thermal X-ray emission concentrated in the SNR rim, with 
filaments that have been taken as evidence for electron 
acceleration^ l4ll28llT3l. 

VERITAS observations performed between 2008 and 2010 
reveal one of the weakest sources yet detected in TeV 
gamma rays, with an integral flux above 1 TeV of 1.87 ± 
0.51 sta t x 10 -13 cm- 2 s -x , or - 0.9% of the Crab Neb- 
ula emission above the same energy JT]- The TeV photon 
spectrum can be described by a power law: dN/dE = 
C(E/3.42TeV)~ r , with T = 
C = (1.55 ±0.43 sta t± 0.47 



1.95 ±0.51 sta t±0.30 sys and 



)xlO" 



L TeV~ 



The source is well-fit by a simple symmetric Gaussian with 
width fixed at the instrument point-spread function (68% 
containment radius of 0.11°) and is thus compatible with a 
point source. The center of the fit, at OOh 25m 27.0s, +64 
10' 50" (J2000), is offset by 0.04° from the center of the 
remnant, with the statistical and systematic uncertainties in 
this position being 0.023° and 0.014° respectively. Figure 
[T]shows that the peak of the gamma-ray emission is slightly 
displaced in the direction of a molecular cloud that is seen 
by integrating over the velocity range between ~ 68kms _1 
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Figure 1: Background-subtracted VERITAS gamma-ray 
count map showing the region around Tychos SNR, 
smoothed with an 0.06 Gaussian kernel. The centroid of 
the VHE gamma-ray emission is indicated with a thick 
black cross HI. Overlaid on the image are X-ray contours 
(thin black lines) from a Chandra ACIS observation! 9 1 an d 
12CO emission (J=l-0) from the FCRAO Survey (gray 
lines) (8). The VERITAS PSF is shown as a white circle. 

and ~ 50kms _1 [17'| and that may be interacting with the 
northeast quadrant of the remnant. While it is tempting to 
interpret this as indicative of hadronic emission associated 
with the cloud itself, not only is the effect not statistically 
significant but the association of the cloud with the remnant 
has been recently disputed[23 1. 

Taken together with observations made at other wave- 
lengths, however, the case for hadronic acceleration within 
the remnant appears strong. Figure [2] shows a pair of sim- 
ple model fits, based on scenarios where the TeV emission 
is dominated by leptonic (inverse-Compton) emission and 
hadronic (pion decay) emission ll22l respectively, and gen- 
erated assuming no influence from the molecular cloud. 
Particle spectra of the form dN/dE = AE""e ( ~ E/Ec) 
have been used, with the spectral index a fixed to a com- 
mon value for both scenarios. The cutoff energy E c 
is permitted to differ as is expected for loss-dominated 
distributions|Q~|. For the leptonic model the derived field 
is ~ 80/iG and the hadronic model has a lower limit of 
~ 230/iG. These values are well above both the usual 
~ 3/xG fields of the interstellar medium and the values ex- 
pected from shock-compression of that field in moderate- 
density environments [6 1 . As a consequence, the VERITAS 
detection of gamma-rays from Tycho can be considered ad- 
ditional evidence for magnetic field amplification with the 
remnant. 

With the models as shown here (fitted only to VERITAS, 
radio, and X-ray data[l|, scenarios where the gamma-ray 
emission is dominated by leptonic or hadronic emission 
are both a reasonable fit, with the preference for hadrons 
resting largely on the fact the field of ~ 80/iG required 
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Figure 2: Emission models fitted to the broadband SED (ra- 
dio, non-thermal X-ray, and VHE gamma-ray emission) of 
Tychos SNR. Upper panel and lower panels show a lepton- 
dominated and hadron-dominated model, respectively. The 
long -dashed curve is the inverse-Compton (IC) component; 
the dashed curve the pion-decay emission. The solid line at 
high energies is the sum of these components while at low 
energies it corresponds to synchrotron emission. 

by the leptonic model is noticeably lower than the values 
of 200 — 300/iG derived from the thinness of the X-ray 
synchrotron emission rim[3|. Along these lines of argu- 
ment, a more recent paper, incorporating both Fermi-LAT 
observations of Tycho and the VERITAS spectrum shown 
here, claims that the only model consistent with the Ty- 
cho broad-band spectrum in this case is one in which the 
Fermi and VERITAS gamma-ray spectrum is produced by 
hadrons lEOlED . It should be noted, however, that regard- 
less of the model used to explain the gamma-ray emission, 
the total relativistic particle energy, which represents a sig- 
nificant fraction of the total kinetic energy of the SNR, re- 
mains dominated by hadrons [ 1 j. 

4 SNR G78.2+2.1 and VER J2019+407 

SNR G78.2+2.1 is a - 1° diameter SNR at ~ 1.7 kpc dis- 
tance. At - 7000 years |[T0| |T6| |T8l it is considerably older 
than Tycho, and is thought to be in an early phase of adi- 
abatic expansion into a medium of fairly low density IT8l . 
Much of the radio and X-ray emission lies in distinct north- 
ern and southern features QUI l25ll . Gosachinskij 1 7 1 also 
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Figure 3: VERITAS gamma-ray image of SNR G78.2+2.1 showing the detection of VER J2019+407. The SNR is delin- 
eated by CGPS 1420MHz radio contours (white) [26]; the solid black ellipse co-incident with VER J2019+407 is the 95% 
confidence ellipse for 1FGL J2020.0+4049; the dashed black circle shows the fitted extent of VER J2019+407. The star 
symbol shows the location of the central Fermi pulsar PSR J2021+4026. The white circle (bottom right corner) indicates 
the 68% containment size of the VERITAS TeV PSF for this analysis. 



identifies a slowly expanding H I shell immediately sur- 
rounding the radio shell which Lozinskaya et at. 0181 sug- 
gest was created by the progenitor stellar wind. A gamma- 
ray pulsar PSR J202 1+4026 is located at the center of the 
remnant. Unlike Tycho, G78.2+2.1's large size makes it 
possible to resolve in gamma-rays. 

The VERITAS Cygnus Region Survey 1 29 1 revealed possi- 
ble emission in the vicinity of G78. 2+2.1, which was con- 
firmed by ~ 20 hours of later follow-up observations. In 
the follow-up data alone, a clear signal with 268 ± 34 net 
counts is detected at a location overlapping the northern 
rim of the remnant, significant at the > 7.5<r level after 
accounting for search trials[2 |. Modeling of the source by 
a two-dimensional Gaussian convolved with the VERITAS 
TeV PSF reveals the source to be extended, with an intrin- 
sic extension of 0.18° ± 0.03° stat ± 0.02% s and a fitted 
centroid position of R.A. 20 /l 20 m 00.0 s , Decl. +40°49'12" 
(J2000)). It should be noted, however, the apparent source 
profile (as illustrated in Figure [3} appears both asymmetric 
and non-Gaussian. 

Figure [3] shows the gamma-ray excess detected by VER- 
ITAS in relationship to the SNR as seen at other wave- 
lengths. The emission seems to follow an arc-shaped struc- 
ture in the radio continuum contours (as seen at 1420 
MHz). The brightest part of the emission is offset by 
- 0.5° from the Fermi LAT pulsar PSR J2021+4026 
(1FGL J202 1.5+4026), but is co-located with the Fermi 
source 1FGL J2020.0+4049 only 9" away. This region of 



gamma-ray emission does coincide with enhanced thermal 
X-ray emission that is suggestive of shocked H I [25 1; cu- 
riously, it lies in a void of CO emission [15|. The TeV 
emission also lies near a region of bright [S II] optical line 
emission within the SNR that is identified as shock heated 
gas based on the [S ll]/Ho! line ratio [19|. 

The relationship of VER J20 19+407 and 
1FGL J2020.0+4049 to G78.2+2.1, and the nature of 
the gamma-ray emission being produced, are not yet 
clear. While VER J20 19+407 could be a PWN related to 
PSR J202 1+4026, the significant angular displacement 
between the two sources, when coupled with the total ab- 
sence of TeV emission at the position of PSR J202 1+4026, 
makes this less likely. Likewise, while it is possible that 
VER J20 19+407 is the PWN of an unknown pulsar in 
the line-of-sight towards SNR G78. 2+2.1, this would 
ascribe the location of VER J2019+407 near the [S II] 
line emission and the enhanced thermal X-ray emission to 
chance superposition. A more straightforward explanation 
of the TeV emission and the features observed in the X-ray, 
optical, and radio continuum would be that the gamma-ray 
emission arises from particle acceleration (either hadronic 
or leptonic) within the SNR shock. High energy electrons 
capable of producing TeV photons via inverse-Compton 
scattering should also produce X-ray synchrotron radiation 
detectable as a non-thermal power-law in the X-ray 
spectrum, for which we currently see no evidence in our 
analysis of the ASCA X-ray spectrum, as shown in Figure 
[4] While molecular clouds are frequently invoked as 
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Figure 4: Top: ASCA X-ray spectrum of the re- 
gion of enhanced X-ray emission coincident with 
VER J2019+407Q. The stepped line shows the fit of a 
non-equilibrium plasma model with kT = 0.5 ± 0.14keV, 
a column density of Nh — (3.7 ± 2.0) x 10 21 and a nor- 
malization of N = 1.8 x 10~ 3 . The line at 1.7 keV is 
identified as due to Si; the fit requires an overabundance by 
a factor of two relative to the solar model. Bottom Panel — 
Residuals from the best-fit model. 

potential sites of hadronic gamma-ray emission, the dearth 
of CO in this region means this is not the case here; the 
target material would have to be H I, for which the H I 
shell surrounding the remnant could provide a plausible 
source. In this scenario, the pre-shock density estimates 
of 1 - 10 cm~ 3 based on the ratio of [S II] A6716 to 
[S II] A6731 line fluxes] 19] appear consistent with those 
extrapolated from the TeV gamma-ray flux. 

5 Conclusions 

Detection of gamma-ray emission from Tycho's SNR and 
from the vicinity of G78.2+2.1 has added to the VERITAS 
catalog a pair of gamma-ray sources with a strong poten- 
tial to shed light on cosmic-ray origin and its relationship 
to supernova remnants. Although the relationship between 
the point-like gamma-ray emission from Tycho's SNR and 
a nearby molecular cloud is not clear, the broad-band SED 
from this remnant favors not only a predominance of ac- 
celerated ions within the remnant but a hadronic origin to 
the gamma-ray emission itself. The extended gamma-ray 
emission from VER J20 19+407, on the other hand, ap- 
pears to track part of G78.2+2.1's shell, which because of 
the remnant's large size is easily resolved. While the total 
lack of CO at this position makes it clear we are not deal- 



ing with a molecular cloud scenario, a hadronic origin to 
the gamma-ray emission from VERJ20 19+407 is not ruled 
out. The lack of evidence for a non-thermal component at 
the position of VER J20 19+407 would favor hadrons, and 
preliminary indications are that the density of H I, possi- 
bly due to interaction with the wall of a cavity blown by 
the progenitor stellar wind, is sufficient to produce the ob- 
served gamma-ray flux in a hadronic emission scenario. 
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